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Section 1 
INTRODUCTION 

Scale model testing which requires simulation of the prototype plume 
shape generally employs a n cold flow" test approach in obtaining the scale 
model plume. Air is most commonly used since air supply systems are 
readily available at most test facilities. However air does not exhibit the 
desired flow properties under all test conditions which then in such cases 
necessitates the use of other simulant gases. Freon 14 (CF^) has been 
found to exhibit expansion properties similar to those of the Space Shuttle 
main engines and thus is a candidate for use as a simulant gas in Space 
Shuttle scale model test programs. 

Previous investigators (Refs. 1 and 2) have found these gases not to 
behave as an ideal gas (constant ratio of specific heats, etc.) nor obey the 
perfect gas law when expanded supersonically. Consequently, analytical 
studies which treat the supersonic expansion of these gases should include 
the effects of the non-standard equation of state as well as the variation of 
specific heat ratio with pressure and temperature. 

Following the work of Refs. 1 and 2 FORTRAN computer codes have 
been developed to calculate the thermodynamic properties of these gases 
for isentropic expansions from given plenum conditions. Thermodynamic 
properties for air are calculated with equations derived from the Beattie- 
Bridgeman nonstandard equation of state (Ref. 1). Thermodynamic properties 
for Freon 14 are calculated with equations derived from the Redlich-Quang 
nonstandard equation of state (Ref. 2). 

Utility of the programs for use in analytical prediction of model rocket 
nozzle flow fields was enhanced by arranging card or tape output of the data 
in a format compatible with a method- of- characteristics (MOC, Ref. 3) com- 
puter program. The nonstandard equation of state programs thus enhance 
the accuracy of the method-of-characteristics solution of supersonic flow 
fields.. 
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Section 2 

DISCUSSION OF THE NONSTANDARD EQUATION 
OF STATE PROGRAM FOR AIR 


The nonstandard equation of state program for air consists of seven 
subroutines which perform the necessary calculations to expand air isen- 
tropically from given plenum conditions to a given pressure. Plenum tem- 
perature and pressure are input to establish an isentrope along which the 
expansion will be made. Along this isentrope the program iterates for tem- 
perature at points located at pressure increments of DELP. At each point 
thermodynamic properties are calculated with equations derived from the 
Beattie-Bridgeman equation of state (Ref. 1). The thermodynamic data is 
arranged into a format compatible with the method-of-characteristics (MOC) 
computer code (Ref. 3) and may be output on cards or a FORTRAN tape unit, 

The input instructions and description of the input FORTRAN symbols 
are covered in the program input section. Interpretation of the output is 
discussed in the program output section. 

2.1 PROGRAM INPUT 

This subsection contains a detailed description of the program input 
as follows; 

• Detailed input guide 

• Detailed description of the input FORTRAN symbols. 
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2.1.1 Detailed Guide for Input Data 
Card No. 1 


F or mat 

.1615 (right adjustecl) 


Column 

Parameter /’ 

Identification 

1-5 

I WRITE 0 

No intermediate printout 


■ ' , — ■ ; , j 

Intermediate printout 

6-10 

ITAB Number of entropy cuts 

11-15 

IUNIT 0 

No output for MOC program 


7 Punched output for MOC program 

N FORTRAN unit for tape output 
for MOC program 

Card No. 2 


Format 

5D10.6 


1-10 

PREF 

Reference pressure corresponding to 
reference entropy and enthalpy, psia 

11-20 

TREF 

Reference temperature corresponding 
to reference entropy and enthalpy, °F 

21-30 

SREF 

Reference entropy, Btu/lbm-°R, 

31-40 

HREF 

Reference enthalpy, Btu/lbm 

41-50 

EMMAX 

Maximum Mach number to be considered 

Card No. 3 ■ ■ 

Format 

5D10.6, 4A6, 111 


1-10 

PTU 

Plenum pressure, psia 

11-20 

.:'gr 

Plenum temperature, °F 

21-30 

DELP 

Pressure increment at which calc ula- 


tions are made, psia 


3 
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Column 

Parameter 

Identification 

31-40 

PMINN 

Minimum expansion pressure, psia 

41-50 

PS 

Plenum pressure for calculation of 
second table, psia 

51-74 

HEADER 

Gas identification for MOC input 

75 

1ST OP 

1 Discontinue calculations with 
immediate case 


0 Continue calculations with 
new plenum conditions 


2.1.2 Input FORTRAN Symbols 


rmbol 


Description 


Units 


IWRITE This flag controls the output of intermediate N/A 

calculations during the iterative cycle for a 
convergent solution. If IWRITE=0, there is 
no intermediate printout. If IWRITE=1 inter- 
mediate printout will be generated. 

ITAB This flag determines the number of entropy N/A 

cuts calculated by the program. If ITAB = 1 
their e will b,e one entropy out. If ITAB =2 
there will be two entropy cuts . If the thermo- 
chemistry data is to be used in a flow field 
containing shocks, ITAB should be set equal 
to 2. 

IUNIT The method of output of the program's thermo- N/A 

chemistry data is controlled by this flag. If 
IUNIT=0, there is no output generated for the 
MOC program. If IUNIT = 7, the output will be 
punched on cards in a format appropriate for 
input to the MOC. If IUNIT= N (N/0, 7), the 
output will be written on FORTRAN tape unit N. 

PREF This variable is the pressure at which the 

reference entropy and enthalpy are evaluated. 

TREF This variable is the temperature at which the 

reference entropy and enthalpy are evaluated. 

SREF A reference entropy for change of ontropy 

calculations. It has a value defined for air 
at a pressure of PREF and a temperature 
of TREF. 


psia 

°F 

Btu/lbm-°R 
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Symbol 

HREF 

EMMAX 

PTU 

TTl 

DELP 

PMIN 

PS 

HEADER 


ISTOP 


Description Units 

A reference enthalpy for change of enthalpy Btu/lbm 

calculations. It has a value defined for air 
at it pressure of PREF and a temperature 
of TREF 

The maximum Mach number which will be N/A 

considered in the calculations of the isen- 
tropic expansion from the given plenum 
conditions. 


The plenum pressure from which air will be psia 
isentropically expanded. 

The plenum temperature from which air will °F 
be isentropically expanded. 

The pressure increment for succeeding calcu- psia 
lations from plenum conditions. 

The minimum static pressure which will be psia 


considered in the calculations of the isen- 
tropic expansion from the given plenum 
conditions . 

Plenum pressure for calculation of a second psia 
table of thermodynamic data. 

Identifies the set of computed gas properties N/A 
when thermochemistry data is punched on 
cards or written on a tape for subsequent 
input to the MOC program. The MOC pro- 

f ram compares column for column , this 
eacler ana the gas identification data input 
to the MOC while searching a tape file for 
the proper set of thermochemistry data. 

This flag is a run control parameter. If N/A 

ISTOP = 1, the program discontinues cal- 
culations with the immediate set of plemum 
conditions. If ISTOP = 0, the program cal- 
culates another isentropic expansion with the 
next set of plenum conditions input via card 
No. 3. ISTOP allows for calculation of isen- 
tropic expansions from numerous plenum 
conditions with one computer run. 
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2.2 PROGRAM OUTPUT 

This subsection contains a description of the output scheme utilized by 
the program and a description of the error messages printed out by the pro- 
gram. 

2.2.1 Description of Program Data Output 

The data output of the nonstandard equation of state program for air is 
organized into three horizontal lines for each pressure increment at which 

4 

calculations are made. Data output includes plenum conditions, local fluid 
and thermodynamic properties and local Mach number and temperature 
assuming air to behave as an ideal gas. 

The following guide consists of numbered flags on a sample printout 
sheet (page 8) which correspond to numbered comments listed below. 


1 

M 

Mach number 

2 

R 

Gas constant for air, ft-lbm/°R-lb-mole 

3 

GAMA 

Ratio of specific heats for air 

4 

T 

Static temperature, °R 

5 

P 

Static pressure, psia 

6 

A/A* 

One -dimensional expansion ratio 

7 

H 

Enthalpy, lbm 

8 

V 

velocity, ft/sec 

9 

PGM 

Mach number assuming perfect gas 

10 

PGT 

Temperature assuming perfect gas, °R 

11 

PC 

Plenum (total) pressure, psia 

12 

TC 

Plenum (total) temperature, °R 


6 
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13 

ENTROPY 

Entropy, Btu/lbm- R 

14 

HC 

Plenum enthalpy, Btu/lbm -°R 

15 

RHO 

Density, lbm/ft^ 

16 

SS 

Local speed of sound, ft/sec 

17 

RHO*U 

•y 

Mass flow per unit area, Ibm/sec-fr 


2.2.2 Description of Program Error Message Output 


1. DUE TO PROGRAM LIMITATIONS, THE PROGRAM IS 
ASSUMING MINIMUM PE= 

This statement occurs when the program is 
unable to solve for temperature at a pressure 
increment. The previous pressure is assumed 
to be the minimum expansion pressure. 

2. NO CONVERGENCE AFTER 100 TRIES 

This statement occurs when the program is 
unable to converge on a solution at a given 
pressure after 100 iterations. The quantity 
generated by iteration 100 is accepted as the 
solution and the run continued, 

3. ITSUB WILL NOT CONVERGE 

This (statement originates in subroutine ITSUB 
which provides iteration control for any function. 
The statement indicates that the program is 
unable to converge on a solution. 
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Section 3 

DISCUSSION OF’ THE NONSTANDARD EQUATION 
OF STATE PROGRAM FOR FREON 

The nonstandard equation of state program for Freon 14 (OF^) consists 
of nine subroutines which perform the necessary calculations to expand Freon 
isentropically from given plenum conditions to a given pressure. Plenum tem- 
perature and pressure are input to establish an isentrope along which the ex- 
pansion will be made. Along this isentrope the program iterates for pressure 
at points located at temperature increments of DELTT. At each point thermo- 
dynamic properties are calculated with equations derived from the Redlich- 
Quang equation of state (Ref. 2). The thermodynamic data is arranged into a 
forrnait compatible with the method-of-characteristics (MOC) computer code 
(Ref. 3) and may be output on cards or a FORTRAN tape unit. 

The input instructions and description of the input FOF,TRAN symbols 
are covered in the program input section. Interpretation of the output is 
discussed in the program output section. 

3.1 PROGRAM INPUT 

This subsection contains a detailed description of the program input 
as follows: 

• Detailed input guide, and 

• Detailed description of the input FORTRAN symbols . 
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3,1.1 Detailed Guide for Input Data 
Card No. 1 


Format 1615 (right adjusted) 

Colu mn Parameter Id entification 

ITAB Number of entropy cuts 

IUNIT 0 No output for MOC program 

7 Punched output for MOC program 

N FORTRAN unit for tape output for MOC 
program 

Card No. 2 


1-5 

6-10 


F ormat 

5E15.8 




1-15 

A2 

Constant for 

nonstandard 

equation of state 

16-30 

A3 

Constant for 

nonstandard 

equation of state 

31-45 

A4 

Constant for 

nonstandard 

equation of state 

46-60 

A5 

Constant for 

nonstandard 

equation of state 

61-75 

A6 

Constant for 

nonstandard 

equation of state 

Card No. 3 

V ' ' j' 


■ 



Format 

5E15.8 




1-15 

i . . . . 

B2 

Constant for 

nonstandard 

equation of state 

16-30 

B3 

Constant for 

1 

nonstandard 

equation of state 

31-45 

B4 ' 

Constant for 

nonstandard 

equation of state 

46-50 

B5 

Constant for 

nonstandard 

equation of state 

61-75 

B6 

Constant for 

nonstandard 

equation of state 
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Card No. 4 


F ormat 

5E15.8 


Column 

Parameter 

■ j< ‘i 

Identification 

1-15 

C2 

Constant for nonstandard equation of state 

16-30 

U C 3 

Constant for nonstandard equation of state 

31-45 

C5 J. 

Constant for nonstandard equation of state 

Card No. 5 



F ormat 

5E15.8 


1-15 

BLIT 

Constant in nonstandard equation of state 

16-30 

AK 

Constant in nonstandard equation of state 

31-45 

ALPHA 

Constant in nonstandard equation of state 


Card No. 6 


F ormat 

5E15.8 


1-15 

AL4 

Constant used in calculation of specific 
heat at zero pressure 

16-30 

BL4 

Constant used in calculation of specific 
heat at zero pressure 

31-45 

CL4 

Constant used in calculation of specific 
heat at zero pressure 

46-60 

DL4 

Constant used in calculation of specific 
heat at zero pressure 
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Card No. 7 

Format 5E15.8 


Column 

Parameter 

Identification 

1-15 

R 

Gas constant for Freon 14 (CF^) 

16-30 
Card No. 8 

AJ 

Constant used in the nonstandard equation 
of state 

1-10 

TREF 

Reference temperature corresponding to 
reference entropy and enthalpy, °F 

11-20 

PREF 

Reference pressure corresponding to 
reference entropy and enthalpy, psia 

21-30 

VREF 

Reference specific volume corresponding 
to reference entropy and enthalpy, ft^/lbm 

31-40 

SREF 

Reference entropy, Bt\x/lbm-°R 

41-50 

HREF 

Reference enthalpy, Btu/lbm 

51-60 

EMMAX 

Maximum Mach number to be considered 

Card No. 9 

1-10 

TO 

Plenum temperature, °F 

1 1-20 

PO 

Plenum pressure, psia 

21-30 

DELTT 

Temperature increment at which calcula- 
tions are made, F 

31-40 

TMIN 

Minimum expansion temperature, °F , 

— - 
i 

o 

PS 

Plenum pressure for calculation of second 
table, psia 

51-74 

HEADER 

Gas identification for MOC input 
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3.1.2 Input FORTRAN Symbols 


Symbol 


Description 


ITAB 


IUNIT 


A2 

A3 

A4 

A5 

A6 

B2 

B3 

B4 


This flag determines the number of entropy 
cuts calculated by the program. If ITAB = 1 
there will be one entropy cut. If ITAB=2 
there will be two entropy cuts. If the thermo- 
chemistry data is to be used in a flow field 
containing shocks, ITAB should be set equal 
to 2. 

The method of output of the program's thermo- 
chemistry data is controlled by this flag. If 
IUNIT=0, there is no output generated for the 
MOC program. If IUNIT=7, the output will be 
punched on cards in a format appropriate for 
input to the MOC. If IUNIT=N (N=0, 7) the out- 
put will be written on FORTRAN tape unit N. 

Constant used in the calculation of the non- 
standard equation of state. A2 = -2.16259 for 
Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

Constant used in the calculation of the non- 
standard equation of state. A3 = 4.404057E-03 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

Constant used in the calculation of the non- 
standard aquation of state. A4 = 1.921072E-04 
for Freer 14 (CF^) in Eq. (Al) of Ret. 4. 

Constant used in the calculation of the non- 
standard equation of state. A5 = -4.48 1049E -06 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

Constant used in the calculation of the non- 
standard equation of state. A6 = 5. 838 8E -07 for 
Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

Constant used in the calculation of the non- 
standard equation of state. B2 = 2.135114E-03 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

Constant used in the calculation of the non- 
standard equation of state. B3 = 1.2828 18E -05 . 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

Constant used in the calculation of the non- 
standard equation of state. B4 = -3.918263E-07 
for Freon 14 (CF^) in Eq, (Al) of Ref. 4. 
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Units 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 
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Symbol 

Description 

Units 

B 5 

Constant used in the calculation of the non- 
standard equation of state. B5 - 9.0623 18E -09 
for Freon 14 (CF^).in Eq. (Al) oi Ref. 4. 

N/A 

B6 .:/ . 

Constant used in’ the calculation of the non- 
standard equation of state. B6 - -9.263923E- 04 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

N/A 

C2 

Constant used in the calculation of the non- 
standard equation of state. C2 = -18.94113 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

N/A 

C3 

Constant us ed in the calculation of the non- 
standard equation of state. C3 = 0.539776 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

N/A 

C5 

Constant used in the calculation of the non- 
standard equation of state. C5 = -4.836678E-05 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

N/A 

BLIT 

This constant is'us'ed in the nonstandard 
equation of state. BLIT = 0.0015 for Freon 
14 (CF 4 ). 

N/A 

AK 

This constant is used in the nonstandard 
equation of state. AK = 0.00976 for Freon 
14 (CF 4 ). 

N/A 

ALPHA 

This constant is used in the nonstandard 
equation of state. ALPHA = 6.61199 for 
Freon 14 (CF^). 

N/A 

AL4 

Constant used in the calculation of specific 
heat at zero pressure. AL4 = 0.0190458084 
for Freon 14 (CF^) in Eq. (A4) of Ref. 4. 

N/A 

BL4 

Constant used in the calculation of specific 
heat at zero pressure. BL4 = 3.00892783E-04 
for Freon 14 (CF^) in Eq. (A4) of Ref. 4. 

N/A 

CL4 

Constant used in the calculation of specific 
heat at zero pressure, CL4 = -1.30237441E--07 
for Freon 14 (CF^) in Eq. (A4) of Ref. 4. 

N/A 

DL4 

Constant used in the calculation of specific 
heat at zero pressure. DL4 = 1 . 96802894E-1 1 
for Freon 14 (CF^) in Eq. (A4) of Ref. 4. 

N/A 

R 

Gas constant used in the nonstandard 
equation of state. R = 0.1219336 for 
Freon 14 (CF 4 ). 

ft-lbm 
°R -lb -mole 
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Symbol 

Description 

Units 

AJ 

This variable is used in calculations of the 
nonstandard equation of state. AJ = 0.185053 
for Freon 14 (CF^) in Eq. (Al) of Ref. 4. 

N/A 

TREF 

This variable is' th 6 temperature at which the 
reference entropy and enthalpy are evaluated. 

°F 

PREF 

This variable is the pressure at which the 
reference entropy and enthalpy are evaluated. 

psia 

VREF 

This variable is the specific volume at which 
the reference entropy and enthalpy are 
evaluated. 

ft^/lbm 

SREF 

This variable is a reference entropy for 
change of entropy calculations. It has a 
value defined for Freon 14 at a temperature 
of TREF, a pressure of PREF, and a specific 
volume of VREF. 

Btu/lbm-°R 

HREF 

This variable is a reference enthalpy for 
change of enthalpy calculations. It has a 
value defined for Freon 14 at a temperature 
of TREF, a pressure of PREF, and a specific 
volume of VREF. 

Btu/lbm 

EMMAX 

The maximum Mach number which will be 
considered in the calculations of the isen- 
tropic expansion from the given plenum 
conditions. 

N/A 


3.2 PROGRAM OUTPUT 

This subsection contains a description of the output scheme utilized 
by the program and a description of the error messages printed out by the 
program. 

3.2.1 Description of Program Data Output 

\ 

The data output of the nonstandard equation of state program for Freon 
is organized into two horizontal lines for each temperature increment at which 
calculations are made. Data output includes plenum conditions and local fluid 
and thermodynamic properties. 
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The following guide consists of numbered flags on a sample printout 
sheet (pp. 18 ) which correspond to numbered comments listed below. 


1 

PO 

Plenum pressure, psia 

2 

TO 

Plenum temperature, °R 

3 

VO 

Specific volume for plenum conditions, 
ft^/lbm 

4 

CVO 

Specific heat at constant volume for 
plenum conditions 

5 

CPO 

Specific heat at constant pressure for 
plenum conditions 

6 

gamma 

Ratio of specific heats for plenum 
conditions 

7 

ENTROPY 

Entropy for plenum conditions, Btu/lbm-°R 

8 

ENTHALPY 

Enthalpy for plenum conditions, Btu/lbm-°R 

9 

SPEED OF SOUND 

Local speed of sound for plenum 
conditions, ft/sec 

10 

TEMPERATURE 

Static temperature, °R 

11 

VOLUME 

Specific volume, ft /lbm 

12 

PRESSURE 

Static pressure, psia 

13 

CSUBV 

Specific heat at constant volume 

14 

CSUBP 

Specific heat at constant pressure 

15 

gamma 

Ratio of specific heats 

16 

RHO 

3 

Density, lbm/ft 

17 

ENTHALPY 

Enthalpy, Btu/lbm 

18 

ENTROPY 

Entropy, Btu/lbm-°R 

19 

SOUND SPEED 

Local speed of sound, ft/sec 

20 

VELOCITY 

Velocity, ft/sec 


16 



LMSC-HREC TM D390169 


21 

MACH NO. 

Mach number 

22 

COMP FACTOR 

Compressibility factor 

23 

RHO*U 

Mass flow per unit area, lbm/ft -sec 


3.2.2 Description of Program Error Message Output 


1. ITSUB WILL NOT CONVERGE 

This statement originates in subroutine ITSUB which 
provides iteration control for any function. The 
statement indicates that the program is unable to 
converge on a solution. 

2. NO CONVERGENCE 

This statement originates in subroutine STATE or 
subroutine SOLVE. In subroutine SOLVE the state- 
ment indicates that the program will not converge 
on a solution for the local temperature. In sub- 
routine STATE the statement indicates that the 
program will not converge. 


17 
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Appendix A 

LISTING OF THE NONSTANDARD EQUATION 
OF STATE PROGRAM FOR AIR 
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GOT 09p 

ITS03300 


" ) 1 IP 

3 3* 


3S 

SAVE 1 1 1 •* 

1 TS03H00 


?;i n 

3 1 • 



SAVE 1 S I -EOF X 

ITSC3S00 


OP i 37 

3S • 



X-a-SavE 1 2 1 

ITS03600 


r .'i 1 3? 

36* 

C 


OF ONE VARIABLE 

1 TS03700 


" o i s 2 

37. 

C 


IF OP* (-FUNCTION Am 1 CM IS DRIVEN TO ZERO 

I T SO 3 e 00 


• : i i J 

ifl* 

C 


(A (•VARIABLE Art I C m IS ITERATIVEL? SOLVED for 

1 TS03900 


niiW 

39* 

c 


• SAVE 1-PROGRAM CONTROL 

I TS01000 

o O 

"'-IT? 

to* 

c 


SAVE 1 2l*X INCREMENT 

t T SOI 200 

£ 


t| • 

c 


SAVE (3)*C0UNTEh DENOTING NTrt ITERATION 

1 T SOI 300 


• ' » \ J 

H 7 • 



Gu I09J 

I TSOttOO 


- in 

")• 


SC 

1 T 1 ME ■ 3 

1 TSOtSOO 

§ £ 

~ •; i » S 

•|t. 



IF IFOF * »GT .C.00IG0T070 


# £■« 

r ' 1 1 ? 

IS* 



IF IF OF AC*. LE. F OF *1 GOT OSS 

I T SOt 7 00 

jO t-rj 

"r i * l 

16* 



SAVE I2)*-| .0;*SAVE 12) 


C 

‘ • i * ? 

17* 



X.X.2.DO*SAVE 1 2) 


s 

"r i * 1 

*1* 



GO T 09" 

1 TSOSOOO 


" ' | t 

IV. 


ss 

SAVE 16) -X 

I T SOS 1 00 


""Its 

S • 



S* VE I 7 1 *EoE X 

1 TS0S200 

"< d 

• 1 * 6 

SI • 



X-A.SA VE 1 2 1 

1 TS0S300 


: I •• 7 

s?* 



G0T09C 

ITSOStOO 


1 ‘•0 

S3* 


70 

1 T IME-t 

I TSOSSOO 


ISI 

St. 



lFIFOF*.LT.0.O0)GOTO7S 



"Cl*! 

ss* 



SA VE I <t >«X 

1 T S OS 8 00 


* > | s*« 

S A • 



SA VE C S ) "FOE X 

1 TSGS900 

> 

mss 

S 7 • 



GoTOR? 

1 TS06P00 

i 

C : 1 '• A 

SR* 

• 

7S 

SAVE (6I"X 

1 T S 36 | 00 

00 

?0 1 s ? 

S 9 • 



SAVEI7I*EoFX 

1 T SO 6 200 


m •; i i ? 

6"* 

C 


HICK NEA GUESS E OR X ACCORDING TO TYPE CALCULATION 

I TS06300 


c > so 

6 1 • 


ec 

XaSAVCMI-SAvEISlMISAVElAltSAVEMII/ISAvEm-SAVCISn) 

1 T 5 0 6 t 00 


••0 1 A 1 

*2* 


9C 

IF I NI.GE. NT I ME SI GOTO 100 

I TS06S00 


. A ) 

S3* 



N | *N I . 1 

I TS06600 


C ! !> *• 

fct* 



GOTO 1 20 

I TSP6800 


" " ; a s 

SS* 


ICC 

1 T 1 ME ■ 6 

I TS06900 


r : l " a 

66* 



GOTO 1 20 

I TS07000 


"" 1 1 ? 

*7* 


1 10 

1 T 1 M* .S 

1 T SO 7 | 00 


» 0 

6 P • 



SA VE ( 9 ) *X 

1 TS07200 


O' l T 1 

69. 



saveisi*eofx 

1 T SQ7 300 


'">1 * ? 

/r* 



SAVE C 6 1 •* 

1 TSOTtOO 


' 0 1 T J 

71* 



SAVE 1 7 ) *EoEX 

I TSOTSOO 


C3I 71 

72* 


120 

CUNT tNUE 



001 7 s 

73* 


SOI 

EONMAT 1 1 X ,SIS) 



0" 1 7 6 

71* 


SC2 

FuRmATI IX,10I2.S«/1X ( 80|2.S) 



POST? 

7 S • 



Y* X 

1 T SO 7 7 00 


03200 

76* 



return 

I T SO 7 8 00 


PO/OI 

77* 



end 

1 T SO 7 9 00 


END OF COMPILATION: no diagnostics. 


LMSC-HREC TM D390169 


if U» , S I T S'JH . I T SUB 

n s o ti* - 02 / 1 v/tr-i j:3i :2s 10,1 


SUBROUTINE ITSUB EnTRY POINT C002II 


STORAGE USED! COUt I I I 000230: 0* T A I 0 ) 000021 1 BLANK C0HM0NI2I 000000 


F*TERNAL PEffRiNCES 

Q 

^ 22 n 00 3 


(BLOCK, NAME) 


NEWR2S 


82 

00CR 

NKRH 3S 











STORAGE 

ASSIGNMENT 

iblock, Type 

, RELATIVE LOCATION, NAME 1 







ID ' 

ss 

000 1 

CC0026 

IOL 

9001 

0C016** 

I09L 0001 300167 

1 IOL 

0001 000177 

I20L 

0001 

0000 A 1 

30L 

£ 0 

?"0 j 

'j99?a2 

3SL 

OCOI 

909071 

SOL 0901 9C0 1 1 3 

SSL 

0001 000122 

70L 

0001 

000136 

7SL 

S H 

0^01 

09? 1 a2 

pol 

cod 

900 1 S3 

99L 0000 0 OCOOOo 

roFx 

0900 0 000002 

FOFXCK 

0000 

000013 

INJPt 


opto 1 

000036 

N 1 

0000 0 

P0C39R 

X 







00 1 

1 p 1 

I* 


SUBROUT 1 NE 

ITSUB 

(FOFY,T,SAVE ,CONV,NTIMES 

.itimei 







OOICI 

2 • 

rn l c 1 

3 • 

9? 1 ? 1 

A* 

, 991 PI 

s» 

vO 99 193 

6 • 

99 I « A 

7* 

rr. irs 

• a 

cr 1 *s 

9a 

" " 1 9 6 

I 0 » 

On | 9 7 

1 1 • 

P j 1 9 7 

1 2 a 

CCIIO 

1 ;• 

90110 

| Ra 

oo 1 1 n 

ISa 

C 9 l 10 

16a 

Mil) 

17a 

90110 

1 P» 

C?l 19 

19a 

- - 1 10 

2 P • 

| 1 2 

2 1 • 

99117 

22* 

"9 1 1 3 

7 3 • 

ruiR 

2 R • 

C'!l IS 

2 *» • 

*- 1(6 

2 6 • 

*01 1 7 

21 a 

' ; 1 2 1 

Vk a 

-0 1 7 3 

2 V • 

nn | ?S 

)"• 


C This SUBROUTINE PR0VI0ES ITERATION control for ant function 
C SAVE (R-7|aST0RAGE LOCATIONS FOR X AND FOFX 

C iCONV ) -CONVERGENCE CRITERIA 

C InT|HESI-HAX NUMBER OF ITERATIONS 

OOUBLERRECISIONFOFT ,r,SAVE,CONV,FOFX,FOFxCK,X 

OIRFKSIONSAVECBI 

F0FxCK-SAvE(P> 

C FOfY ANO T ARE DUMMY INPUT ARGUMENTS 

FOFX-FOFY 
XaY 

C CHECK f OR CONVERGENCE 

IFIAHS (FOFX l-CONv«LE«0* )GOTO| J(J 
C (TIME CONTROLS THE TYPE CALCULATION TO BE PERFORMED 

C I T|ME»l .FIRST TIME THR0U6M 

C ITIME>2,P0S FIRST TIME THM0U6H 

C I T I ME a 3 ,NEG FIRST TIME THROUGH. 

C ITlMEaR, SOLUTION IS BRACKETEO 

C ITIME-S, SOLUTION HAS CONVERGED 

C ITlME«c>, SOLUTION A I L L NOT CONVERGE 

GOTO! »0, 3o,Sp,7e » , I T I ME 
c initialize 

ic n 1 • 1 

I T 1 HE a ? • 

F Of XCK«f UFA 
SAVE I 8 1 aFOFXCK 
If IFOFA.LT.O. IGOT oSO 
30 If (F0FX.LT.9. IGOTO70 

If If Of »C<.Gr .FOFX » GOTO 3S 
SAVr I 2 ) .»SAVE I 21 


ITS00200 
I TS00300 
I TSOOmOO 
itsoosoo 

I T SOOfcOO 
I TS00800 
I TS00900 
I T SO I 000 
I T SO I |00 
I TSOI 200 
I TSOI 300 
I TSO I 500 
I TSOliOO 
I TSO I 700 
1TS9I BOO 
I TSOI 900 
ITS02000 
|TS02|00 
I TS02200 
I TS02300 
I TS02H00 
ITS02S00 
1 TS02600 
I TS02700 
I TS02B00 
I TS02900 
I TS03O00 
I TS03 | 00 


LMSC-HREC TM D390169 



00 I 7 6 

3|* 



X>A-2,*SA vE 121 

1 T SO 3 200 


n Qt 27 

»• 



C.OT090 

1 T SO 3 300 


| jo 

3 3 • 


3 5 

5A ve 1 S 1 -X 

1 TS03S00 


1 1 1 

3** • 



SAVF ISI-FoFX 

I T SO 3500 


•■' 1 3 2 

35* 



XaX-SAVE (21 

1 TS03&00 


r ~ i j 2 

36* 

C 


OF ONE VAKlAHLE 

1 TS03700 


mw 

3F • 

C 


IFOF* »-FunCT|0N AmICM IS DRIVEN TO ZERO 

I TS03800 


"’u? 

38 • 

C 


IXI-VAKIABLE «M|Ch is ITERATIVELY solved for 

ITS03900 


r-’ii? 

3 V • 

C 


1 SAVE 1 -PROGRAM CONTROL 

I TSOSOOO 

S 2 

Tl M 

r <•• i ■» 7 

so* 

C 


5 * v t I 2 ) • X INCREMENT 

1 T50S200 

T*J 

. H 

'“'M? 

•* 1 • 

C 


SAvL 1 3 ) -COUnTEk DENOTING NTH ITERATION 

ITS0S300 

% o 

rti'i 

92* 



60T090 

I TSOSsOO 

§2 

~ 1 !«• 

S3* 


5C 

ITIME-3 

1 TSOSSOO 


**' 1 iS 

s s • 



IF IFOF *>&T • 0 • 1 GO T y F 0 

I TSOsfcOO 

~ & 

r - 0 1 3 F 

*3« 



IF (FOFACX.LE.FOFX IG0T05S . 

I T SO s 7 00 

c 'v 

oc 1 * 1 

s 6 • 



SAVEI2I--I.*SAVE(2I 

I TS2N800 


"**1*2 

S F • 



X.A.2, • 5 A vE 1 2) 

1 TS0S900 

F ° 

r : i s 3 

S 8 • 



GO T 090 

1 TSOSOOO 

b 

r > i <4 s 

S 9 • 


55 

SAVEI6I-X 

1 TSOSlOO 

S© 

C .0 1 * S 

t*» 



SAVEIFI-FOFX 

1 TS05200 

Ml'* 

SI • 



x-a»savE 1 2> 

1 T SOS JOO 


**" 1 s 7 

S 2* 



GOToVo 

1 TSOSsOO 


00 1 SO 

S3* 


70 

ITIME-S 

I TSOSSOO 


00 1 *• i 

SS» 



IF (FOFX.LT.O. IG0T07S 

I T SOS 7 00 

> 

, S 3 

ss* 



SAVE IS l«X 

1 TSOSSOO 

i 

00l«.** 

S6* 



SAVE 1 5 1 »F OF X 

1 T 505900 

►— 

00ISS 

S F • 



GO (000 

I T 5 06000 

o 

C'llS* 

S 8 • 

• 

75 

SAVE < 6 1 *X 

1 TS06 | 00 


''"IS7 

S9» 



SAVE(7»*FoFX 

1 TS06200 


!*•■' 1 S7 

60* 

c 


PICK NE* GUESS FOR X ACCORDING TO TYPE CALCULATION 

1 T S 0 6 3 00 


nu 1 ) 

6 1 • 


80 

XaSAVEIII-SAVEISIslISAVElAI-SAVElRH/ISAvE 1 71-SAVE CSI 1 1 

ITS06S00 


: -» * h l 

*2* 


90 

IF INI .GE .NT I ME SI GOT 01 00 

ITSOASOO 


00 1 6 7 

*7* 



N | • N | ♦ | 

1 TSOAiOO 


00 1 *«•* 

as. 



GO TO 1 20 

ITS36U00 


00 1 6 5 

*S* 


ICO 

1 TlME-6 

1 TS06900 


00 * A* 

66* 



GO 1 0 1 20 

1 TSC7000 


00 1 *7 

*7* 


1 1C 

ITIME-5 

1 T SO 7 | 00 


r. : \ | TO 

6 H • 



SAVEIStaX 

1 T 50 7 200 


001 FI 

69* 



SA»E I 5 1 »F OF X 

ITS07300 


fill 

/C* 



SAvr.i6i*x 

1 TS07S00 


00 1 7 J 

71* 



SAVE 1 7 1 *F OF X 

ITS07S00 


COI 7R 

72* 


120 

CONTINUE 



"01 FS 

73* 



T • X 

ITS07700 


Out 7* 

7 S • 



RETURN 

1 T SO 7 BOO 


0.0 1 F F 

7S* 



EnO 

1 T SO 7 9 00 


E N o OF COMPILATION; NO DIAGNOSTICS. 


LMSC-HREC TM D390169 


Vfgl ,s ITIRaT.ITERAT 

«•»« I I A -32/ I V/ 7*-l 3 I 31 5 2a 10, » 


SUBROUTINE HERAT EnTRT POINT 000101 


storage used: codeiii ccoiso: dataioi oooq7o: blank comhoni2I oooooo 


(ohkoi blocks: 

0503 INPUT 000030 

caternal references (Block, name: 


OQOn 

l tsub 

ooos 

Nf H<* 2t 

MO* 

NAUUS 

09U 7 

n*02* 

09 10 

NtKHJI 


STORAGE ASSIGNMENT (BLOCK, TYPE, RELATIVE LOCATION, NAME I 


CON 1 


cooa 1 3 

1 CL 

0001 


000071 

isl 

0001 

0000*5 

20L 

0000 


000027 

2*F 

o 

o 

o 

0 

000012 

A 

5303 

0 

C03CC* 

AO 

0003 

0 

'*930 1 N 

a 

009 3 

0 OOOOlp 

BO 

0003 

0 

0000' * 

C 

0000 

0 

00002 A 

C ON V 

Muo 

u 

9 0 90 2 2 

For 

0090 

1 

00902* 

IBR 

0090 

0000*9 

1NJPS 

0003 

0 

00002 A 

ocv 

0003 

0 

000022 

OBAN 

030C 

0 


pa 

9003 

D 

99590* 

PE 

9093 

0 OOOOOO 

PTl 

0003 

0 

000020 

R 

0000 

0 

000002 

SAVE 

0C03 

0 

9PC926 

TE 

0993 

0 

050092 

TTI 













00 1 o I 

i • 


SUBROUTINE HERAT I PR,T,Er*E2,E3,RH0,l TIME) 

0 0193 

2 • 


COMMON/ INPUT /PTl , TT| , PE, AO, BO,*. , B , C , R . 0G AM , OC V , TE 

9310* 

3 • 


OOUBLEPRECISIOnPH.T.EI ,E 2,E3.RM0,PA,SAVE,F0F,C0 NV,«,PT| ,tti ,pe, 

"9 1 r.t 

*» • 


• AU ,nn , a ,B ,C ,OGAM , oc V ,TE 

53 1 ?S 

"> • 


0|MENSI0NSAVE(8I 

9C | r * 

*• 


RhU* p R / >' R • T | 

O’l 1 9 7 

7 • 


1 BN* 1 

MHO 

A • 


SAVE ( 2 1 • • i 09 

9 i. 1 | 1 

V* 


CUNV-.OC5O50C9IOC 

00112 

10* 

10 

Pa*RHO*R*T*(I , 09*HH0* IE l*HM0*(E2*RH0*E3l | l 

03 1 1 3 

1 1* 


FoF«P«-PA 

"Cl 1* 

1 2 • 


call i t sum if of, mho, save ,conv , too, ibr ) 

331 1 S 

13* 


GOTO! 10, 10. 13, 10, IS, 20) , IBR 

9911* 

1 **• 

20 

AM 1 TF t*,2SIPR,M,T ,RHO ,FOF 

991 2S 

IS* 

2 S 

Format I 3i]ttH|TSUB *)LL NOT CON Vf RGE ,P , R • T , R«0 ,F0F« ,50 1 2. S | 

Ml ?* 

1 6 • 


1 T 1 MF ■ | CO 1 

“ >l 77 

1 7* 

IS 

Conti nuE 

' f » o 

IP* 


RETURN 

-r i » i 

1 *♦ 


EMU 


f no or compilation; no diagnostics. 


t- 

s 

CO 

n 

I 

X 

H 

M 

O 

H 

s 

d 

U> 

vD 

O 


O' 

NO 



.rO«,S PtN'i»'|Pl«(>»S 

-.J I I A -02/I7/7H-1 J: 31 530 «C,I 


SUHRQUlINf PERGAS EnTRT POINT C03073 


storage usrc: code (i i oootos; dataioi cooojmi blank commonc2> ogoooo 

conhcn blocks: 



> 

I 

» — • 
IS) 


0003 INPUT 0C3030 

EXTERNAL REFERENCES (BLOCK, NAME I 
000m xpuo 

OOOS DSU * T 

OOOA NEHK3S 


S T OR ALL ASSIGNMENT (BLOCK, TYPE, RELATIVE LOCATION, NAME) 


0001 


000060 

IOL 

0001 


sooois 

20L 

0003 

D 

OOCOI 2 

A 

0003 

0 

00000* 

AO 

0003 

0 

COCO 1 0 

BO 

0003 

0 

ecooi* 

C 

COOO 

0 

no"o06 

FACTI 

0000 

0 

000010 

F ACT 2 

0003 

0 

000022 

GAMA 

0000 

0 

OOOOOR 

Gm 1 06 

0030 

D 

ooeooo 

GM 1 02 

0000 

0 

000002 

GM 1 02 | 

oro 3 

0 

OOCC2R 

OC V 

0003 

0 

90000“ 

p 

0003 

0 

COOOOo 

PO 

0003 

0 

000020 

R 

0003 

0 

CC 0002 

TO 

cooc 

0 

3000 1 2 

XM 










rotoi 

i* 

00 1 C 3 

2* 

00103 

3 • 

00 1 OM 

N* 

OOICS 

S* 

Oil"* 

6* 

| 0 2 

7 • 

«0 1 1 3 

8 • 

roll* 

»• 

0 2 1 1 2 

IC« 

f'C 1 1 3 

1 1 • 

''Oil “ 

12* 

001 IS 

13* 

00116 

I M* 

Or | 20 

IS* 

OSI 21 

16* 

0.122 

17* 

1 23 

IB* 

"31 2M 

1 »• 

9 1 2S 

20* 

C | 2 6 

2 1 • 

r - | ft 

22* 

r c ( )0 

2 3* 


subroutine pergasipm.pt, ierror.pmim 

double precision pm ,pt ,po , to.p , ao ,bo , a ,b ,c ,r,gama ,oc v ,t ,gm 102 , 
l6Hl021»GM|0GiFACT|,FACT2tXM,FACT3 

C 0 M M 0 N / I N P U T / P 0 , T 0 , P , A 0 i B 0 , A , B , c • R • G A M A • OC V • T 

Gm | 0 2 •(GAMA-l.l/2. 

Gm 10 21 * I • /GM I 02 
GhlOG - (GAMA. | , I /GAMA 
20 CONTINUE 

FaCTI • I (PO/PI**GM|OGI-l. 

F A C T 7 - FaCTI*GMI02I 
IM «USGRT (FACT2) 

F AC T 3 • ( I »*GM|02*XM**2> 

T • TO*l I./FACT3) 

IF IT.GT.70.0D9IG0 TO 10 
P • P * C.07IDQ 
ILHROR • I 
PM IN • P 
GO TO IV 
10 CunTInuE 
PM • X M 
PT ■ T 
RE 1 URN 
END 


0003 

0 

00001 “ 

B 

0000 

0 

0000 1 A 

f ACT 3 

0000 


000026 

INJPS 

0003 

D 

000026 

T 


LMSC-HREC TM D390169 


F NO OF COMPILATION: 


no DIAGNOSTICS 





LMSC-HREC TM D390169 



■ro*,s sInm.sinh 

MSO II* - 02 /l»/ 7 R-lltli:U 10,1 

function SINK ENTRY POINT 000030 

storage used: coociii ccocm oataioi 0000211 blank connonizi 000000 


EXTERNAL REFERENCES IBLOCK, NAMtl 

OOOJ OF X P 
OOOh nEHRJ* 

STORAGE ASSIGNMENT (BLOCX, TYPE, RELATIVE LOCATION, NAME I 

0000 0 OOOOC 3 El COOO 0 OOOOOS E 2 0000 OOCOU INjPf 0000 R OOCOQO SINN 0000 0 000001 SNMAftC 


> 


OCIOl 

1 • 

00103 

2 • 

09 1 OR 

3* 

90I0S 

R» 

c? 1 "6 

s« 

9C|ol 

*• 

90 1 1 9 

7« 

001 1 1 

*• 

ENO OF 


function sinhia) 

DOUBLE PHECISION A , S NH A R6 , E | , E 2 
SNHARG ■ A 
E | >OEXP(AI 
E 2 aOEXPI-A) 

SINN • C.S* IE I "E2 I 

RETURN 

END 


compilation: no diagnostics. 


LMSC-HREC TM D390169 


.* o» ,S |CA. , |C*t 

"Su I 1 A -02/ I *»/7**-l 3 : 3 I : 33 (0,) 


subroutine ical entry point coilSl 


s t o h a u t used: coueui ooiissi OATAIO) 0003*»7| blank C0MM0NC2I 000000 

common blocks: 

oooj input oncio 

upon f acton oncu 

0?0S S A C T OOOOO** 

EXTERNAL REFERENCES (BLOCK, NAME) 

000* |TE«aT 

Oro; APDI 

00 1C OS 1 NH 

00 t I NAOUS 

0" I 2 N I 0 I % 

00 I 3 N I 0 2 » 

oo i ** oe*p 

OOIS KPOO 

00)6 NlRRJi 

> 

STORAGE ASSIGNMENT I BLOCK , TYPE, RELATIVE LOCATION, NAME I 

LP 


0001 


ccooss 

1166 

OOCI 


000213 

1336 

0000 


000 1 **p 

|S»F 

0001 


000** **2 

1 7 36 

0001 


000377 

SOL 

0001 


fnft; | 

*•01. 

OOPO 


00CIR2 

sor 

0000 


OOP 1 7* 

60F 

0001 


OOOSOi 

6SL 

OOOQ 


000210 

70F 

000 1 


001077 

*CL 

0001 


001103 

ROL 

0003 

0 

000012 

A 

0003 

0 

000006 

AO 

0003 

0 

00001*1 

B 

0003 

U 

OCCOIO 

60 

C003 

0 

000016 

c 

ooos 

0 

OOPCOO 

CAPE 

0000 

0 

oeooi** 

CAPEN2 

0000 

0 

000016 

CAPE02 

POOS 

0 

PPPCO 2 

CAPF 

0000 

0 

00007** 

CAPF | 

0000 

0 

000076 

CAPF2 

0000 

0 

000100 

CAPF 3 

0000 

0 

0QCC62 

CAPS 

pr>00 

H 

0001 37 

C APGN2 

rcec 

0 

00C06C 

CAP602 

0000 

0 

OCOOSo 

CP 

0003 

0 

ooooso 

CPFACI 

0000 

0 

0000S2 

CPFAC2 

PPOO 

0 

pr?"*** 

Cv 

ooco 

0 


CWIB 

0000 

0 

P0??3a 

EA 

0000 

0 

000036 

EB 

0000 

0 

O000**0 

EC 

PPOO 

0 

ro 

CN20EN 

0000 

0 

ocooo** 

EN2NUM 

0000 

0 

00"0 l 2 

E020EN 

0000 

0 

000010 

E02NUM 

0000 

0 

000020 

El 

0000 

0 

C"0022 

E? 

cooo 

0 

"0002** 

E3 

000? 

0 

oe"i 2 e 

G 

ooco 

0 

000026 

GF AC 

0000 

1 

00013** 

S 

0000 

1 

0001 36 

1 COUNT 

ooco 


0003 10 

I NJP* 

000.0 

0 

0C00S6 

I PROP 

000** 

0 

000000 

1 1 

0000 

1 

0001 3S 

J 

0000 

0 

••OOOaR 

JF AC 1 

nor, n 

0 

"00066 

JF AC 2 

0000 

0 

00007*3 

JF AC 3 

0000 

0 

030072 

JPROP 

000** 

0 

000006 

JT 

OCOR 

u 

POCOC 2 

Jl 

0000 

0 

OCOI' * 

KPR0P 

COON 

0 

OOCOIo 

kt 

000** 

0 

00000** 

«l 

0003 

0 

00002*1 

OC V 

0"03 

0 

00002 2 

0 Cj A M 

0003 

D 

O0"00 • 

«*E 

COO 3 

0 

OOOOOn 

PT 1 

0003 

D 

000020 

R 

0000 

0 

000030 

RMO 

"000 

0 

?r*c i ? ** 

RT 1 

eou3 

0 

000026 

TE 

0003 

0 

000002 

TTl 

0000 

0 

000032 

T 1 

0000 

*> 

000000 

ACP6N2 

0000 

0 

COPCO 2 

XCPG02 

0000 

0 

OOOOR2 

2 














"cj«o «o I Agnus t I c • the name capgni appe«rs in a dimension or type statement but is never referenced. 
*>0101 I* SUBROUTINE ICAL I I AH I TE .RHOI .gamma , I TIME > 

“~I»J *• UOUrtl E PRECISION XC PGN 2 , K C PGP 2 , EN 2NUM , E N20E N , E 02NUH , EO*OE N , 

'i i 3» * E APt N2.C APF 02 

i UUOHLE I'RECISIun E I ,E2 .CA.GF AC .RHU.T I ,EA ,Eb |EC ,2 ,CVV 18 ,CV,CPFAC I . 

I CPF A n, CP, I PROP, Capon | ,CAPb02,C«PG.JFA'l , JF A C 2 , JF AC 3 , JPROP , C A PE . 
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<■-106 

6 


* l ' s 

7 


":i"S 

8 


1 ‘ 4 

9 


l - 7 

1C 


. . ; i o 

1 1 


•"MO 

1 7 


- : i 0 

1 3 


1 • 3 

1 - 

%% 

'•Ml 

* l • 7 

IS 
1 6 


. t 1 

1 T 

& 

: i s 

1 M 

: • * 

1 9 

A ^ 

' " 1 7 1 
■ 1 ’1 

2 1 
2? 

l 

23 


r i /s 

2 s 

* l 7 A 

? \ 


' 1 .*7 

26 


• ! » 7 

27 


1 TO 

2 M 


' - ( ,.r 

29 


- i M 

3’ 


1 17 

31 


•" 1 ts 

32 

> 

1 

H-* 

' . ' 6 

"■•I H 

33 
3 - 

r > , <r 

3 S 

O' 

' • i - 1 

36 


' t-7 

37 



38 


- 1 «. S 

39 


" ! t - S 

SP 


" ’i 

s 1 


"l“S 

s 7 


1 6 

S3 


, to 

ss 


" 1 . 1 

SS 


. , - | 

S 6 


l •- 7 

S 7 


r i 

s A 


"1*3 

sv 


• - | s 

SN 


■ • : r s 

SI 


1 S6 

S 2 


• -* | N 7 

S3 


• 1 

Ss 


• 6 1 

SS 


: >.7 

Sfc 


• 7 

S 7 


. • 3 

Gf 


l '• - 

SV 


' ! • S 

6C 


1 

6 | 


2CAPFI,CAPF2,CAPF3,CAPF,KPR0P.RT|.G.PTI,Tti,PC.A0.i* 0,A,B.C,R,0GAM. 
30C W , Tt a I I , Jl f K | , JT ,> I ,RmO| .GAMMA 
COHMON/ | NPUT/PT | .TTt.PE.AO.eO.A.P.C.R.OGAM.OCV.TE 
Common, f AC TO"/ I I » Jl , * I . JT ,KT 
COPMON/SFACT/CAPt.CAPF 
DIMENSION l<Tt|»)i(il6) 
c calculation of dens:ty 

C UNITS LF VARIAoLES t ARE F 1 • • 3 /L B • MOL E 

C UMTS OF VARIABLES G A"E in--2/lbf 

t: • HO- 1 AO/ 'W»TE I l-C/TE-O 
I S I FOHMATI |A.70I?*SI 

£ 2 • A<>« A/ ( R • TE ) -PO-B-BO-C/ TE-» ) 
f 3 • I Ft 0 • B • C I / T C • • 3 
Do 10 1-1.6 
X T I ( X I ■ R • T E • • 1 
10 CON T t Nut 

(.1 I »• -t 1/HTI I I I 

GI2I- I2.C*F|-*?-l2)/RTI (2) 

6(31- (S.f'-IE?-EI«-2l-E3)/RTI (31 

o < s > • ( 6 .P*F I -t 3-3. 9-E2--2- ( 7.g«C I--2I -({.O-E I ••2>3«0*C2 * ) /HT M H I 
G ( » I ■ T.t»ir?-t3-‘f.0-E|*lE2--2*EI-E3l*l6.0-EI*-3l-l2«0-E2-EI--2l » 

\ /HT||S| 

GI6I- I 6 « C • I I IO*C*C|--2I-I3.S-E2--2»2«0-EI-E3I«(||.0-E|--*I) 

* • I 2 • 9*E I ••2-S.0-E2l-2.0*E2--2» -E3-I H . 0 -E 3-7 3 • 0-E I-E2I I /RT| 1*1 
gf'«: • i.p 

O' 20 1-1,6 

GF AC > GF*C * Gill -PE »• | 

20 CGNTINUC 

c the real gas densiit is nor computed 

CALLITEFUTIPE.TE.EI .E2.E3.RM0l ,|T|ME» 

RhO-HHOI 
Tf-TE 
EA-E I 
Eu-E 2 
EC -E 3 

c Tut ratio of specific heats calculation 
c vihrational contribution for air 
IF I { Time . £Q, II GO TO 30 
2-.C68S6100 

:wlrt » I I 706.850C/TE 1/SlNMC 30Ss.900/TE » >--2 

• ♦ I C?MS.2?D0/TEI/SINH|20R6. |D0/Tei»--2 

Cv • OCW ♦ CvwIH*(*»0-2-C-RM0/TE--3I • I 1 , O-BO-RHO/2 • O-BO-B -RM0-*2 
»/3.C» 

CPF AC | -I I I •9-2.0-C-RM0/TE--3I-I I »0-BO*RhO-BO-B-RHO-- 2| I • • 2 
C PF A C 2 • I I • 8* ?• C*E 1 •RH0-3.0»E2-RMU--2*R»e-E3-RH0--3l 
CP - C V ♦ 2- ICPFAC I /CPF AC2 I 
GAMMA • CP/CV 
GO TO SC 
30 CONTINUE 

Gamma ■ ogam 
so CONTINUE 

c Calculation of iprop 

1 1* « 0 p ■ <GAmma/0Gam)-I|.0«2«0-C| -RhO* 3.C-E2-RH0--2-R «0-E 3-HH0--3 I 
It l I -R I TE • l*T , o I R R | TE 16, SOI I T IMF ,E I ,E2,E3,<GfUI , J ■ I ,6) , GF A C ,RH0, 
SGAMMA , | PROP 

lFII>R|Tt.GT.M.ANi>.IT|HE.GT.IlRRlTE(6,63)CVVlB,CVfCP 
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M02M3 

AC- 


•••.;? I 1 

6 I* 


r 2 1 ? 

6R • 


-i2l? 

6S- 


y ■> ? 1 3 

A6- 


' •/ 1 R 

6 7. 


• 2 1 5 

6A- 



4»« 


" ?! ? 

70 • 


' ' :o 

71 • 


r • . 2 1 

;?• 


** 2 2 1 

■ I 

7 3- 
7R- 

sf 

. 2 

< 2 1 

■ 

,'S 

»S- 

76- 

77* 
7 8- 

D rO 

2 6 

7 9- 

fZp- 

• >.-27 

BO* 

'fr o 

• .• *C 

8 ; • 


. > 1 

«?• 


"-.M2 

83- 

•'“2 32 

HR* 


"12 0 

MS- 


"02 

86- 


0 -• 7 x S 

« 7- 


r ; 1 6 

8 H * 


-'2 T7 

89- 

> 

2-0 

Vtj. 

1 

o i '» i 

9 1 • 

►— 

: i s i 

92- 


"•3/K | 

93- 


' r ? 

9 R • 


"02S 3 

9S« 


- S3 

-*6- 


•'MSI 

97- 


" “ 2 *• 6 

9 0- 


r ,«.< s 7 

99« 


r "/60 

ioc- 


*07A 1 

»ai- 


? 2a i 

1 02 • 


Po.- il 

103- 


" 2 A 5 

IOR- 


f-.'AA 

105- 


- - 2 A 7 

104- 


"f.2 70 

1 37- 


7? 1 

1 OR • 


r w ;? 

109- 


r ■ . » J 

1 IT 


0"/ 7 R 

III* 


" 2/5 

1 1 / • 


* . : ; 6 

1 13- 


""277 

1 1 R* 


<:•••• calculation of jprop factor using capg 
Icoont • e 
6 s continue 

C • • • Capg FOR aIR FROM REFERENCE EQ. 78 

C a P G N 2 • <-R770.9n0/T| I/I I.C-DE(P<6I09.70Q/TI I I 

CAPG07 ■ ( -8S7 .290Q/T I I / I I .P-DE *P I 9092. |0«/T I I I 

CAPO ■ C APGN?*c*PG0? 

JIACI » RhO- (MO- t 2.0* AO/ I R-T I I |.R.O-C/T I ••] } 

JF AC? - R m u • • ? • I ( 3.P-AO-A/I2.0-R-T | | l-S,o»ttO«C/(2,0»T|..3)-B0»BI 
JIAO • »HO« »3- ( 2. Q-bn-B-C/T | • . j i 

JPnuP - I • C •» I COGAR-I.0l/OGAHI»(CA * G ♦ JF A C | ♦ jF A C 2 ♦ JF A C 3 I 
C---» CALCULATION OF KPROP FACTOR USING CAPE And capf 
C-* CAPt F NOH REFERENCE CS. 7S 
* C P 6 N 2 ■ C*PGN2/0 .780B8O.0 
ALPf.02 • C APGO2/0.2CVS00 
EN2NUM • 0t«P(*CPGN?| 

EN2PEN ■ I .O-OE AP | -6 109. 700/T I | 

E0 2’4UM • 0£ API *CPG02 I 
E J20FN ■ I . R-OEXP | -HR92. IOO/T | I 
C A PE N 2 • |tN2NUN/tN2DEN I •• I c.78"8a00l 
C A PE 0 2 ■ iEU?NJP/F 02DEN I •• CP. 2095001 
Cape • Capcn?*capeo? 

C«» CAPF FkOM RFrCHfNCl EU. .38 

CAPFI • Kmu* ( HO-2 . ?»C/T I "*3 I 

CAPF? - KH0--2-I (UO-C/T|»-3I-BO.B/2.0> 

CaPFJ ■ RmO* • 3* 2 a" -bO-B -C / I 3 <0* T I • • 3 I 
CAPF -UEAPlCAPF | *CAPF2-C APF3I 

A P R OP - ICAPE/CAPF )•( I .o*ea*Rho. E 8*RH0**2«EC*RH0«*3I 
lF<|6H|Tt.GT.0»*N|TE(6, 70 J RMO , C A PG , JPR0F , C A PE , C 7.PF , K PROP 
5C FORMAT ( |H’, 37 m|N SUBROUTINE ICAL F OR ITERATION NO. ,I3,9ME(t-3I • 

I • J02P. I 2 t / I H | 9 HG I 1-61 • , 6020,12, /IH ,6 h GF 1< C • , 020 . I 2 , SH RM0*,02(J 
2.12, 7h GAhma-, D 2C. 12, 7H I PROP- , r>21. I 2 I 
60 Format I |m?,7m C V V I b- , 020 • I 2 , 2* , AM CV-,02p.!2.?X,RH CP-, 020.121 
70 FORMAT! Imm.Sm RMO- , D 20 . I 7 • 6M C A PG- , 0 2C . i ? , 7H JPHOf** »020» I 2 , 6M CAFE 
I ■ , U2G a I 2 , 6H CAPF- ,D2f. I 2./7H KPROP- ,020* I 2 I 
IF « I T I ME ,C«. I . ANO. ICOUNT. EQ. I I &0 TO 80 
|| - | PROP 

J| • JPROP 
K| ■ KPROP 

lF||T|ME.Gr.n GO TO 90 
RMO • P T I / I R-TT I I 
ri • TT I 

Ea • BO- I AO/ |R-T I I l-C/T | --3 
Em ■ AO-A/ I R-T I I -B0-B-80-C/T | e-3 
EC • BO-II-C/T I - • 3 
ICUUNT-I 
GO TO 6S 
60 C ON T I NUf 
JT • JPROP 
KT * KPROP 
90 CONTINUE 
Rt I OP N 
F NO 


cro of compilation: i diagnostics. 
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Appendix B 

LISTING OF THE NONSTANDARD EQUATION 
OF STATE PROGRAM FOR FREON 14 (CF 4 ) 


- . » >. 1 n 

..a - ?/:■',/ 7i - 7 i : ?t : it ic.) 

»<!', PH*J 5 ‘-AM ' ... 

S?..*‘Ot. oSlOl (<.Uf ill CC 1 C 17 1 CATAI&) COSlbc! 0LANX COmmoM2> 032C33 


I ^ U *• i I LHti-tl.b iLLIlf, I.AMf.) 




-Oi 

5* *TL 


y«. 

^ t 5 

r « 

i-Ci 

L j l* 


. . 

l .5 7 

S'wl.U 
l > »>AL 

oS 

^ S' 

i .. IL 

S( LVL 


; i ) 

HiS 


. -i / 
. i 3 

I TiUB 
.• J I.TH» 


l:li 

fc •« L O A 

|D *tJ 

• 1 b 

h t l . 2 % 

£ £1 

. 1 fr 

l.liJJ 


• wl 7 

hi « T J 


• .-’I 

5LA T 


. / 1 

M>H2» 

; 7? 

r..-t u* 


. . 2 J 

1. 1 o 1 » 

. . 

< _7i 


** 


y 


i .25 I.S10P6 


aSSIOWKLUT ( OlOO > ‘ T YU , RELATIVE LOCATION, NAME) 


OH 


C C H 7 7 6 

1 CCA 

CLCC 


ucsccc 

1 C 1 E 

0301 


000S15 

1 1 5 L 

0001 


030610 

1 30l 

□ 001 


000531 

1 3 1 L 

3 ..01 


cc:t*b5 

1 3i L 

CCCiO 


CCS 1 1 1 

l 1CF 

COCO 


CCS J 15 

USE 

OoCC 


035 1 20 

1 RUE 

OOC1 


003617 

1 6 l L 

LwO I 


C ( L <>06 

1S7L 

CCOI 


C C C 6 6 6 

l&’L 

0001 


0001 7 s 

19L 

OwOl 


OUOM 7 3 

1 9UL 

0001 


000675 

1 9 1 L 

l >Uti 


t.Lbl C2 

i’CLf 

cccc 


ccs i : i 

20lf 

coco 


OC 50 <2 

2 0 2 (■' 

0300 


OC S0 1 3 

2y3f. 

U030 


0051 10 

20**E 

i wCl* 


C (.5 107 

/ct r 

CCIC 


ccsicc 

20*F 

505 1 


033773 

2 | C«- 

G 330 


00S122 

2 1 bp 

oOOl 


000272 

262b 

• vOI 


( ( C 173 

3L 

ccco 


C C 1 7 7 b 

3c<-r 

coco 


CoS 1 05 

JC 1 E 

Olu 1 


03 1 0 1 3 

303l 

UOw 1 


COO 7 1 0 

100L 

L ..0 i 


rc C7t6 

1 1 CL 

ICGl 


CCtb77 

M 1 7 & 

0031 


0CC4S6 

HSlb 

ULC 1 


C3LM5 7 

bL 

uooj 


C30163 

7 L 

1 ot 


ILtl 12 

VCL 

(■CC 7 

R 

CCCC2S 

A J 

00 J 2 

R 

03C0I6 

ak 

0302 

R 

C3C01 7 

AlMmA 

U0U2 

K 

003020 

AL1 

.10 

► 

tLl7i6 

AM 

CCC2 

R 

cccc:c 

A 2 

30 32 

R 

OOC 3 7 1 

A3 

OuG2 

R 

3 3001)2 

A *4 

0002 

R 

OCOOO 3 

A 5 

(•2 

h 

L( l COM 

a6 

CU 7 

R 

CCCC |1> 

blit 

CO 32 

R 

3CC02 1 

BL 1 ) 

3w0 2 

R 

coooos 

h2 

0002 

H 

033006 

b 3 

•. . P7 

h 

os. toe? 

1 *1 

CCC.2 

f 

icrcu 

c s 

0032 

R 

03051 1 

B 6 

OuuU 

K 

00*«7h1 

cf 

0002 

H 

000022 

CLl 

ct.J’f 

K 

U 1771 

cor.v 

C C t->C 

b 

OC M72C 

CONVR 

COCO 

H 

O \j *i 7 6 0 

Cp 

OwvO 

R 

03i7jb 

CPO 

0000 

R 

031757 

CV 

v?:CC 

H 

U.173M 

eve 

uu 

b 

c : c : i * 

C 2 

0032 

N 

COCCI 3 

C 3 

OlL 2 

H 

035011 

C 5 

OCCO 

R 

001726 

OCLTT 

1 _r/ 

h 

LI 3023 

IL1 

L( :b 

b 

• * r • » • 
VCvbv% 

C jR - - 

033 *3 

R 

uocc 1 c 

E M 

Cvi-G 

H 

001 723 

L M M> A X 

0007 

M 

3U0C00 

LTHAL 

1 _w( : 

h 

Cl i 773 

f U 

Cl ct 

b 

C C M 7 A 1 

6 AM 

33 w 3 

R 

3Ci 7 36 

Ga«3 

G^cD 

R 

kj j 1 7 6 0 

faAMZ 

0000 

R 

001761 

pan 

« ,ul 

it 

Cl *t 7 1 b 

c. A M 

cccc 

b 

c;-7si 

OR 

39»«* 

R 

001762 

M 

JwvO 

R 

031772 

HC 

ocoo 

n 

03171 1 

HEADER 

l I L C: 

M 

( l i 7 •• 1 

> C 

1 1 1 2 

h 

CCCC 32 

h*ff 

00 53 

I 

Oil 73 1 

1 

3wc 0 

I 

031771 

I e* 

U0u3 

I 

0(ji76 7 

1 OF 

: i> 

! 

LI i < ->* 

1 b 

Cl C l 

l 

c:*>72 1 

1 TAB 

0 333 

1 

Owl 7 ,? 2 

lyNIT 

Out O 

I 

001 7 d 6 

3 

OOcO 

1 

031753 

L 

i «• 

1 

. r •* 7b2 

t. 

cccc 

A 

C C i 7sfe 

R 

CO-., 

K 

Of i7 >5 

»'0 

JC.O 

« 

33373L- 

MOT 

U0L2 

K 

030027 

PmEF 

► 

' 1 -WO 

1 A 1 s 

lire 

A 

; c*>7 j. 

► 5 

CCwj 

h 

g.li 1 1 

»*%> A 

-A - .4 

K 

0 <3 i2i 

H 

oQwJ 

k 

C3&30J 

Sa ft 


r 

-7 4. J 

•- 1 ► 

i L 

n 

c;i; ,i 

SlOC 

• c * ^ 

* 

. . ti >< 

1 3 * 

. * - 6 

•I 

O-ivob 

syoi , j 

.’ll 

H 

Ww J J J 1 

SKLf 


h 

: l • ‘ i c 

•s 

cccc 

k 

CCS7S1 

T 

Cl»u 

A 

*..17 75 

T L rJ 

•J w w J 


031/27 

l .i l u 

UC30 

>1 

0 U 1 7 2 1 

TO 


» 

• / 7*i 1 

•t « 

till 

► 

; s 7 H 7 

TP 

- -2 

k 

. : 7 6 

T* f ► 

j . t. 

if 

f. .• 1 7 6 ’* 

J 

0030 

If 

301755 

V 
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CO b rc*’733 VC. 


tccc <■ c . 9 7 s ; vp 


0002 « OOOC30 VrEF 


COCO 9 030/79 2EK 



f : »l-i 

> 


r .. i r i 

; 


: t e'- 

2 


er 1C? 

* 


C 7 J £ *. 

< 


Cl I c 

6 


' : i c ~ 

7 


: 1 1 3 

t 


Cell** 

V 


cr ! 73 

1C 


£.179 

1 1 


LL 1 32 

1a 


'*:,I‘.L 

1 3 


C c I ** !< 

1“ 


CLISJ 

IS 


CIS/ 

1 /• 


6 7 

I 1 


■ : 7i. 

1 « 


-* rcc* 

I v 


•- ..-Cl 

;t 


* . r r j 

2 1 


7C-9 

2/ 

» 

; • - t 

* • »/• 

/ .» 

I 

- .'Co 

is 

CsJ 

' ’ ’ C 7 



: ’ic 

26 


' r,7i i 

27 


!>_7I2 

2c 


f*r ? 1 3 

2 V 


t r 2 1 v 

3 C 


7 m 

3 1 


C .23 1 

32 


C C 7 3 2 

33 


c r 2 3 J 

31 


r„7 as 

3b 


CC23S 

36 


C .- 7 3 A. 

37 


C 0 2 3 / 

3b 


Cf 791 

3 V 


ft 2 s 1 

*•0 


rr 7ci- 

*1 I 


CL2S2 

92 


t » s 2 

93 


t 7 *< 2 

99 


'•.7*3 

9 S 


r . *•*. 

9/. 


* • V, A 

* r 


* »l # 

A » 



V 


r 

1 


C 0 » r t r. #i,*3i***.P-5,46,B2.a3|b'«.aS,B6.C2 l C3,C5 

COP rir |.LJ t ( A&. ALPHA. AL9.*.L9 .Cl-F.0L9.N.AJ 
COP rtf- 1 F L F ,F't"EF.V«e F fSREF,rtKEr 

Ul l '•*. It ' SfVVUO lER(SOC) •ZEMl&aci ,GA«Zt500> ,T0T|S00).P0TtS3C» 
• » h i a I r f i •• i 

ctr vp» 2;» i ?‘i*77u«2 • 

ke 1 1 « * ,?t i r 1 1 At ( J unit 
3C0 ECU t I II t 1 C* I 

R t / 1 < ' , 1 r c IA/.aj, A 9 .A 5 .A 6 
ico id r * i iif is.t i 

Hut .i r eJt2i9j..j** t os,a« 

M«.L «•- . lCClC2.C3.CS . . . - — 

M » l i i , I C C * t L I T , a*. , ALPHA ©C 

HE , 1 1*. . i : £ ) *(.*• .uLH ,cL*t ,i>L9 W *2 


r»Ml«L,l?C)K,AJ ptfS 

Ht/ L to* I.* I UP E F .FR£F aVREP .SREF»N*er .EMM AX Q g 
»d cuMini Q ta- 

9 HC/UL.jt I»IC, PC. 0 EI.TT, THIN, PS, mEaDE* fc* 

I’jl FdrFlUllt.il p> ^ 

If lit •L1 »»$!LGTc 3C3 CJ s» 

j-i £ Q 

is-L C E# 

call state i pc. to. voi g ** 

Cali r e at s t to . vo »cvo »cpo . saho i .. ... . . 

St.i t *L* A I TC , vO ) 

SS«St t r-o t T 0 , * C »& AMOI 
PL-I r I F *ETmAL ( To. vo> 

19 Cor.Tir.L-E 
PPM T 2Ct 

SLLc • SCSh . . ................ 

PH 1 1 1 2tt,PO.TO,VO,CVO,CPO.aAMO,SOSR»HO,SS 
PSF a- PC-199, 

Cf -P 0- VC/ t F- TO I 
OASI*IS9E.2«32.l7M»cF/8e»OI 
AP-l .C 

. - .... 

2EF- ( I r • G A SH 
ti A F A I I I » f- A h 0 
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12 .' 

2 3* 

C-.I 23 

2m* 

C » ! 7 3 • 

• 2 S • 

00! 2N 

26* 

C C t 2 *• 

2T 

CC! JV 

2t>* 

CC I 2 <• 

29* 

cr 1 2 ? 

3C * 

Col 30 

31* 

CM 31 

32* 

Cel }i 

33* 

CO 1 3 3 

3*<* 

CC»3‘« 

3S • 

CC ! 30 

36* 

Cl’l 3o 

37* 

CC! 37 

3B* 


j I ,C/<2.0*0! *•?>*■• •0/ , 2»0*02**2! ) 

T0L..,.C.<*^c,.r».C d UT/<i.a-U|.-s)-BLIT/«H.C«0 2 ..m- 

Tr,».i» ■ TCa*T 0 (!*TDC*TCD*T« 3 c*TDF • 

TEA »t.T«*L 06 « 02 /i'j»*l*tB 2 -AK»C 4 !»eKT»*«l« 0 / 0 1 -l.ij^D 2 l 

TEl«lMh3-*<«CJ»t<1>»l0.S'Dl ,# 2-0.S/02**2> 

TEC*T »fil* < > • :/o I • »3-> •0/‘J2*«3 » / J»0 

TEL,*1*U-S-a**CS*LiCT>*»I .C/0 »••*<»» •0/02»«H/h*0 

Hf»T , 66 , tE* v -‘ tAVf<,/Al ' Prl * 

Tl«A». , *»TlA*TE;»*TLC*TEWoT£fl 

ETt'Al»TA*TiJ*TC-Tl**T£ 

KETOM. 

END 


tm, o» compilation. 


NO DIAGNOSTICS. 


w 

I 

vO 
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r » n n 


iSp.tTfc 

•«<,{ I | A -l»2/7U^7*t-iI »C*> 

rU.CTlCI. IS* _ INTM POINT CCC2bt _ 

SToRAGt I'btCj COOtjl) CCC 26 S ; OATaiO) 0000331 dLANK CO*HON( 2 ) 00 q 033 


f 7 TLRNAL UFEPl.f.CIS IOlOCA, n A h C ) 

( CO J EXP 
ICO** ALO& 

COOS MMO* 


SToRaCi ASSIGNMENT (BLOCK, yTPE, RELATIVE LOCATION, NAME) 


w 

I 


LOO 2 

h 

CLC025 

A J 

C t C2 

R 

CC00I6 

AK 

0032 

R 

330017 

alpha 

0002 

R 

030020 

AL 4 * 

0002 


000000 

A2 

LCC2 


ttccoi 

A3 

CCC2 


cccc:2 

A <4 

0332 


300303 

A5 

0002 


03000** 

AO 

0002 

R 

03301b 

0 L 1 T 

1002 

H 

CC002I 

IL*» • 

CCG2 

R 

c.c:;s 

0 2 

0332 

R 

00003* 

03 

0002 

R 

030007 

BM 

0002 

R 

030013 

bS 

t.i.02 

h 

CCLOl » 

l .6 

CCC2 

K 

3C0C22 

CL 4 * 

C032 

R 

030312 

C 2 

0302 

R 

030013. 

-C3 — 

.0002 

1 

.0030 1 *4 

cs 

L o02 

H 

CrtC23 

l L** 

rccc 

■■ 

f r*f * * 

05R 

GO 3 3 

R 

030302 

01 

0300 

R 

030003 

02 

UOGO 

R 

00030 1 

E*T 

CLC2 


LLCCJ2 

».RIF 

tccc 


COCO 1 3 

JNJPS 

3032 


30332/ 

PREF 

0302 

H 

00002*) 

R 

0002 

R 

000331 

s«ei 

f.LOC 

H 

CCCOL’A 

T CMP 

CCC2 

K 

CCCC26 

TREF 

0332 

R 

303033 

VrEF 










r *• » r, 1 

CCIC3 
C.IO 1 * 
C'. It'S* 
L,.JLA 
cr. I L 7 
C I I u 

cr. 1 1 1 
r t i 1 1 

C_l 12 
2^11/ 
C; ! I 3 

• 0113 
'..11** 
: i i *» 
~l I !» 
r. 1 1 A 
f r 1 1 7 


3 • 
««• 
i* 
6- 
7* 
8 • 
9 • 
1C- 
1 1 - 
12 - 
13- 
1 1- 
lb- 
16- 
1 7- 
1 U- 


FUNCTION osr<t*v> 

COMMON A2,A3,A‘«,AS,A6,B2ia3,B < l,dS,B6<C2,C3 f CS * 

COMMON pL I T , Alt , ALPHA » AL9 iilL** »CLH »DL«* ,«• AJ — .. - 

COffCN TNEF,PNCF»VReE*5NEF,mREP 

EaT«L7.P(-AK*TI - « 

0 I »VH F-t-L I T 

02»V*bL I T * ......... ... 

TEMfAL‘i-AL0G<T/TREF>-BL ,, MT-Trtt.E>*CL**MT**2-TR£r**2)'2»0-0L‘»* 

I <T-»3-TMF--3)/3*C 

TE»F-UM’*Aj*AA-lLKT-EXP«-AK»TRE F »»-«C2/02-c3/<2*0*02»*2) 

|♦C^/< , ^.0 , C^•• , ♦) )-Aj*N*ALOg«u2/01 I 

Ul » •11P.P-A j* I J»2-AK*C2*E<T )•« I O/tta-I *0/D| * 

1-Aj«U3-IA»C3-EaT1*<1»0/iJ2-»2-1.0/01--2)» , 2»0 . 

UMP"lENP-Aj*(U5-AA.Cb*EAT)-CI.U / 02-- , «-l.O/ol* , 9)/*»*0-Aj*89* 

1 < I »t./l 2 • • 3- 1 •3/Ol**3)/J*3 

OSF«TtNp-AJ*<96*ltXP» ALPMA»v 1-EXP1 ALPMA-V R EF ) >/AlPHA-SNEF 

PE TOM. . 

CM 


c» lompjl*»;om 


NO OfAGNjbTtC** 
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M C* 


• 3 SOLVL.SULVt 

•'•jt. u* -C2//0>'7*»-2| :^6:b* «c.» 


S».pl»0‘M pil SOL Vf EliTRY POINT cccno 


STyk^Ct OSLO; COL’L(I) CCCJJOi CaT«<C) CCCONJj BLANK cOMMON< 2) 000033 


IaTERNaL » EftHLNCtS (BLOCK* IlAhtl 


0003 


USH 
















COON 


I T S UH 
















COOS 


A 1 NX 







• . . 




.. ... 





L.06 


CAP 
















000 7 


1*1 NN 2k 


.. 

.. 

.1 

* . - «. .. 

- M M 




___ 






COO 


f.l«H3* 
















STu» AOt 

A*-*S | GNltLNT 

(BLOCK, TYPE* RELATIVE LOCATION, 

NAME) 

— ... - — - 

.. - ... 



- ... 

. ........ 

• • . - . . . 





CCOI 


CL 0077 

1L 

. ... CCCC 


c::ci* 

1SF 

0031 

_ OCOOJ* 

*L 

COOL 


_ 03007?. 

3 l 

0002 _ 

00002S 

A J _ 

LC02 

P 

Cl CPU 

AK 

CCC2 


crrci 7 

ALPHA 

0S32 

000020 

Al*« 

0002 


000300 

A 2 

0302 

00300 1 

A3 

LJ.02 


CCC002 

A*» 

CC 02 


o c c c ;• 3 

AS 

• 0002 

00003*1 

A A 

0002 

R 

0003IS 

B L 1 T - ■ 

0002 

300021 

8 L 41 

L L 0 2 


C( cccs 

02 

CCC2 


c:rc:» 

63 

0002 

OwOON 7 

ON 

0032 


000310 

BS 

00C2 

00 JO 1 1 

66 

Gi,C 2 


CCo022 

CL 4 * 

ccco 

K 

&::: 12 

C ON V 

CO 32 

0C001 2 

C 2 

0002 


03001 3 

C3 

00u2 

0000 1 *• 

cs 

L J02 


croc 2 3 

GL 4 * 

ff“C3 

K 

3 ::;cc 

CSR 

0300 

R 30031 1 

01 

OOGO 

R 

0J0310 

CKT 

UP02 

030032 

href 

l (ICC 

1 

Of 1.0 j s 

Ibh 

coco 


OCC 52 7 

1 MJP* 

0032 

- 030027 

prC* 

- 0o02 

— 

-0-JU<)2‘* 

R 

O 0 C 0 R 

u J 0 O 30 

Saul 

LL 00 

N 

CLOG 1 3 

SL 

Cf 00 

R 

03CC 1 <* 

SOF* 

0002 

030331 

SrEF 

0002 


0 3002 A 

TKtF 

0 002 R 

000033 

VREF 


0~ I C 1 
Oc I 03 
3c I C’’ 
C.-IOS 
!C* 
r ic7 

r-c- i i »■ 

r - 1 i » 
C-v I I 2 


1 • 
2* 
3« 

S* 
6* 
7 • 
H* 
V* 


SUbfcOGT | N£ SOLVE<T,y*TP,VP,tiP**oSK> 

CCPPCl. A7.A3,A<«,AS l AA,B2i43,BS l BS,06tcri 

COM. Cl* ellT l A*tAL’’HA«AL < *><3L‘ltCL < l'0l l *'*»l 

COPNOt. tPEt.P«EF.V*CF.SKtF,rlREF 

D|p| I.SION SAVE I 8 I 

V-VM t lp/1 1 •• I I .i / I fiP-| *01 | 

E A 3 »l. AP ( • A K • T I - — 

Oi«v*tr-eLU 

SAVU I !• I • 


. 1 1 J 

10 * 



SAIL «2>*.CI 

r : 1 1 “ 

1 1 • 



COI V • 1 .F -03 

r. i it 

12 * 


2 

sifctSF n ,v 1 

Ojl 1 6 

13* 

- • • 


SOF A«SCSA-SL • 

c c : 1 7 

It* 



CALL 1 (SUE ISCFX , V.SAvE ,CuNV , 1331 

t- 1 2G 

l!>* 



It F-SAVE 1 1 1 

f:rn 

1 A* 



00 I0l2,2(2*2*l(3l*| BN 

V ' 1 22 

» 7* 


J 

CCMIl-LE 

"• 123 

1 *.• 


IS 

FOFiniUCISP NO CONVERGENCE 1 

; • • ; - 

- IN* 

- 

1 

CU T Il.tL 

• • - r 

. • 4 

i r • 



Hit t t. 

'Ml 

i 1 • 



ll.L 
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r<o ci cohpiim ion: no 


w 


DIAGNOSTICS. 
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t -rci*»s PRtl> • PHES 

•’S'- 1 1 A -U2/70/7R-2I :2A5Sr 1C.) 

FUNCTION f-HES ENTRY POINT COCUS 


SToPaGL U31G5 COOLIll corilll OaTaIO) OOOOISt BLANK cOmnoN<2| 00^033 


EATthHAt fcLFtRfcf.CES (bLOCK, NANEI 


W 

OJ 


LUO 3 


E*P 

















... 





CiUOR 


ME RK3» 


















STORAGE 

assignment 

( BLOCK t TYPE 

. RELATIVE 

LOCATION, 

NAME! 





— - 

■ 


, . - 



•-’302 


LC0025 

A J 

CP02 

R 

CCCCI A 

AK 

0002 

R 

00001 7 

ALPHA 

0002 


000020 

AL" 

0002 

n 

003000 

AZ 

t J02 

X 

LLLOO I 

A 3 

CC 0 2 

N 

ccocc* 

AR 

0032 

R 

ooaco 3 

AS 

0002 

R 

00000" 

A 6 

00w2 

R 

03331* 

ML 1 T 

uL C2 


0LCC21 

PL** 

0C02 

h 

CCCCiS 

b 2 

0C02 

it 

cocao 6 

B3 

0002 

II 

0 JOOO 7 

ai 

00U2 

a 

0J30I0 

MS 

0002 

h 

CC-POI 1 

b a 

CCO? 


CCCC22 

CL" 

0002 

R 

COOOI 2 

C 2 

0002 

R 

000013 

C 3 

0002 

a 

sooom 

cs 

cuC2 


0LLL23 

M.N 

LC CO 

P 

CvCCrCl 

IKT. 

0002 

• • 

003032 

HREF— 

OoOO 

— 

000 0(j3 

|*JP* 

0002 


00O027 

PREF 

CLOU 

h 

CLtCOC 

PRES 

0102 

R 

0CCC2" 

R 

C002 


CQC031 

SaEF 

0000 

R 

000002 

TO 

0002 


00302A 

TrEF 

CC02 


CC LO 30 

VhtF 





•• 




— 










ooioi 

4* 


00103 

2 • 

. . 

CCIOR 

3* 


Cu j cs 

R« 


OUI Ob 

%• 


c j. r>7 

■ - A* 

— 

S3 1 1 U 

7* 


OH 10 

«• 


0-3 1 | 1 

*• 


'1 I 1 7 

10* 



Function presit.vi 

CONNLN A2>A3iAl|/<S»AA»S2>b3tB < t«bS ( B6tC2(C3|CS 

COMON bin ,AK,ALPHAfAU‘«iBL 4 t.CL*».Ol.‘t»R,AJ 

COi.t-UN T*EF ,PRCF * VREF t SRtF ,hREF - - - 

CaT»| AP ( •* AK • T I 

T 0 »v-hL! T •- - 

Phi S»h*T/TO*« A2*t2*T*C2*t*T >/TO»«2M A3 *b3*T»C3*EKTI/TO**3M A"a8«» 
IMI/Tt^-HMAS^B^T^CS^EKTI/TO^htlAA^BA^TI^EiPlALPMA*'') 

RETURN 

ENL . 


ENl> 0» COMPILATION! 


NO DIAGNOSTICS. 
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noon *) or o ci o o 


• f OH , i SOU'.l* . iUllNO 

II* -02/20/7*1-21 :2N{S? (0.) 


rtficrioN sound Entry point 0C0I27 


STuKAOE UiLDI CoDLU) 0001**0I DAtAlGl C0002*ll BLANK C0 nM0N<2) 000033 


E»tE«N*L "CFERtllCLS (OLOC», NAPE) 

0303 LAP 
0_0-< SUNT 
0003 hCnh J» 


STOHAbt ASSIGNMENT ( BLOCK • TYPE, RELATIVE LOCATION, N*HE> 



0302 CC002S 

AJ 

C C 02 

r 

C0C01* 

AK 

0002 

R 

0000 1 7 

AlPMA 

0002 

000020 

AL*» 

0002 

R 

000000 

A 2 


LCC2 h LUCCOl 

AO 

C C 02. 

H 

CC0C02 

A1 

0CD2 

N 

00 3CC 3 

AS 

0U02 R 

03000** 

A* . 

3002 

R 

0300 1 S 

BLIT 


0002 CC002 1 

l»L** 

CCC2 

ft 

CtCC;S 

f»2 

0002 

R 

ccooo 6 

B3 

0302 R 

03033 7 

BH 

0002 

R 

030010 

US 


ww02 H LCLLU 

t* 

C CC 2 


00CC22 

C L*l - - • 

C302 

R 

003312 

• C2 -- 

- 0032 - ft 

000313 

Ci 

-0002 

R 

000 J 1 1 

cs . 

Cj 

1 

0302 OCCP23 

CL** 

croc 

R 

tccri i 

EKT 

0CC2 


CC0C32 

HrEF 

0000 

030UI I 

INUP* • 

0003 

H 

03000 3 

PP V 

1 

h-* 

CJC2 CCCC27 

PNLF 

CCL2 

K 

00002** 

ft 

LCwO 

ft 

C03C00 

sound — 

0002 — 

0 0033 1 

skef. 

0000 

H 

0300C2 

To 


0002 CC 00 2 A 

0 

TPI F 

0CC2 


CCCC3C 

VREF 





... 


-■ 



- 



onoi i* 

■ • -■ 

function 

SOUND! T t V 

• gamma ) 

***“ * 

— 








— «• • — — 



U 1 03 
OIO* 
01 OS 
310* 
3 I C 7 
Cl 10 
0 110 
CIlO 
"III 
3112 
n i3 


2* 
3» 
*•• 
S* 
6 • 
7* 
B • 
?• 
IC« 
t I • 
I 2* 


COM. Of. A2.AJ,A‘l,Ab,A6,P2.e3,B < 1,0S l B6.C2,C3,C5 
COMtCI. bt 1 T ,AK , ALPHA . ALS i3LM ,CLN ,Di.H .R ,*J 
COM.cn TREF.PHF .VReF.SREF.mREF 

EKT-t XP |-AK*T ) 

TD-V-ILIT 

ppV»-K*T/T0»»2-2*5*<A2tB2«T»C2*tKT)/T0**3-3.O»(A3*B3*T»C3»EKTI/ 
lTD*» < <-‘*.C»lA*t*B‘l*T)/TD««i i -S.C»lA5*BS , T*C5»EKTl/T0*«B*ALPHA* 
2lAt«b**|l»UI ULFma»V) 

SOLI. I -SUNT l-0ANNA»PPV*v»*2>«*B. 0*232** 

NETLM. . . . . _ 

tr.L 


r.Hy Of COMPILATION: NO DIAGNOSTICS. 
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rnn,:» i t st;*: . 1 Tsuu 

1 t A -U2/X0/7N-2 l : 2 / : C J 1C.) 


SUBROUTINE I T Son tMM KIM C- C C 2 *• «» 


s)oa*ct usi.o: cudkiij etc-?*-*: lataici coco2c; blank c°mmoe*<2| oooooo 


MTEKNaI KLEEREuciS JhLl'L) . l.M’-t) 

OOCS *.FRN2* 

C0G9 NF.kHS* 


SI0«A6l ASilCNl'LHT lOLOC'. TTPt. ILLATIVE. LOCATION. NAPE )_ 


0301 

0003H7 

IOL 

CL01 

CC.C2 1 2 

I COL 

0C01 

ooczts 

1 1 OL- 

0001 - - 000225 1 20L 

0001 — 

00006 1 

sol 

CUOI 

COO 1 00 

SSL 

1 »0 1 

ecru 7 

SOL 

COOI 

0001 30 

SSL 

OOri 000I3 7 7 0 L 

3001 

000162 

7SL 

0001 

QC3 1 66 

DOL 

COOI 

CCS 1 76 

70L 

coco 

R C00002 

FoEX 

OCOO R 00000 1 FOEXCX- 

— 0000 

- 000010 

! N j. * 

COO'J 

I 00300*1 

I TIME 

LCOO 

i ecccc: 

Til 

' CCOJ 

R COCCOJ 

X 






C3I0I 

W : ' i •: i 

c - 1 C i 

u» : : i c i 

: t o i 
C3 1 0) 
Clio** 
OC I Ob 

co ios 

CjlO* 
02102 
C0I07 
C2I 10 
CO I l 2 
C01 12 
OCll? 
CO I I 2 
COI I 2 
Col 12 
CO I | 2 
CO) 12 
Co I I J 

oo 1 1 j 
oci |1 
tins 
:;ii* 

1 1 7 

• 

. ! J ^ 

. I - 


I* 

2* 

s* 

*i • 
S* 
6* 
7 • 

»• 
9 • 
10* 
1 1 • 
I ’• 

1 3 • 
11* 
IB* 
16* 
17* 
Ik- 
17* 
2C* 
21* 
22 * 

2 3* 

2H« 

2S* 
26 * 
27* 
2d • 
If 
JJ» 




• SUltkOl-TlUC 1 T SUU (F0FY f Y*SAVE,C0NV,NTlMCS) ... 

this suepouth-e e-rgvides iteration controi eor *ny function 

SAVE I *1-7 < -STChAGt LOCATIONS fOR a AND FOf* 

ICMIV >*CONVf K6Er.CE CRITERIA 

If.TiM Sl-E-AX I.UNLER OF ITERATIONS • 

U I N|.)< S I ONS A V E ( 6 I 

Ml “SAVE I 3 i ♦. I 

Foe XCKabAvE I B ) 

FO r Y AN(, Y Akt &UHHY ; NPUT ARGUMENTS-* 

FOF X ■! OF Y 

ChLCK V OR COI.VEK.fHCE 

lElAbS )FOFX l-Col-V.Lt.O* IGOTCIIO 

mnC-SAVC) I »♦. I 

- ITIHL CONTROLS The type CalCULATIONMO-BC P£RfORNET> 

IT1HI *| »F IRST TIME *H*0U6M 

ITIML-Z.I'OS E I K S T T | HE TMROUOH — — 

|T|lil*3.l<t(i FIRST Tint Through 

I T 1 IiF* • SlitOT | CM | S fcKACFt'EO - 

I T I H L • S 0 L U T 101. hAS CONVERGED 

ITICEaR.SOLUTlU; * I LL DOT COI.vEROE 

1*01 lit 1C 3u.SC*7C» , 1 TIHE 

Initialize - 

10 M|* | 

I T | |tl *2 ... 

foe aci *r or x 

SAVl.lolal GFXCK .... 

]| l|M | X.l 1 .C • ) LCK.SC 
33 IE l| M X *| I ICC i 07C 


iTsooioa 

ITS00200 

j T5C33Q0 

ITSOOROO 

ITSCOSOO 

itso3:,0o 

ITS00700 

I TSUOBOo 

ITSC0700 

ITSS1000 

I TSOI 100 : 

1 TSC12O0 

ITSC1300 

ITSOIHOO 

ITSCI&OO i 

I TSl I BOO 

1 T SC l 7 C3 

I TSOI 800 

| TSw I 900 

1 T SO 2030 

ITSC2100 

t TSC22J0 

ITSC2330 

ITSC2N09 
ITSD2SCO 
I T SC 26u0 

I TSCa/CQ 

| I Su2B J J 

ITSC2YJ0 


Hi M aCi .ML *E OE / K.0T03S 


I TSC3L-.C 
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'.H7i 

31 • 



5»VC«2l*-l .*SAVf • 2 1 

1 T SC 3 1 30 



c-u/ 

32* 



A*A-2.*SAVL 1 r 1 

ITS&J2C3 



:;i jo 

33* 



sot ovu 

.. . ITSC3300 



com 

31* 


35 

S A vC ( t ) • X 

I T SC 3 lOO 



C2132 

3!»* 



SAVLlSl*! Of X ... .... 

_ . 1 TSG3S0C 



C2t 33 

36* 



X-a-SAVI. C 2 1 

I T So 3 A Oj) 



to 133 

37* 

C 


OF Of.E V * k 1 AbLE 

2 T SO 3 7 00 



CC» 33 

3b* 

C 


(Fill A)-fIU<CT ICS I'HlCH IS ORlVEu to zero 

1 TSC3adO 



c r 1 33 

3v* 

c 


IAJ-VaRjalll SmICh is ITe«aT1Vi.L v SOLVED for 

. . ITSC39 jQ 



C.3I 33 

sc* 

c 


<SaVC >-F*OGRAM COUTkOL 

I T SC 9303 



••Cl 33 

S|* 

c 


SaVL ( 1 I at T IKE • . 

1 TSC110Q 



V. 133 

•*2* 

c 


S A V ( ( 2 ) ■ A 1NCI-I hE*'T 

I T Sa 1200 



C3J33 

*3* 

c 


SAVEI3»*C0Ur.TLk DENOTING I.Th 1TEKM10N .. . 

ITS313U3 



:o : 3*i 

HH* 



GOTOVC 

IT SOU g3 



C3I3S 

IS* 


SO 

» T I I.E -*3 

ITSCiSlia 



CCI 3* 

••A* 



Ir'IM X • GT • 0 • 1 GOT 070 

I T Sul 600 



CrilaO 

17* 



II iiui xr*r.(.r .fof x i f.oTc*.s s _ 

I TSOI 7uO 



CsU2 

ib* 
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